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Abstract 

A microscopic analysis of the non-dissipative force Fnd acting on a line vor- 
tex in a type-II superconductor at T = is given. We first examine the 
Berry phase induced in the true superconducting ground state by movement 
of the vortex and show how this phase introduces a Wess-Zumino term in the 
hydrodynamic action Shyd of the superconducting condensate. Appropriate 
variation of Shyd gives F^^^ and variation of the Wess-Zumino term is seen to 
contribute the Magnus (lift) force of classical hydrodynamics to F^d- Because 
our analysis is based on the true superconducting ground state, we are able to 
confirm and strengthen earlier work by Ao and Thouless which examined the 
Berry phase arising in an ansatz for the many-body ground state. We also 
determine Fnd through a microscopic derivation of the continuity equation 
for the condensate linear momentum. This equation yields the acceleration 
equation for the superflow and shows that the vortex acts as a sink for the con- 
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densate linear momentum. The rate at which momentum is lost to the vortex 
determines the non-dissipative force F^d and the result obtained agrees iden- 
tically with the Berry phase calculation. The Magnus force contribution to 
Fnd is seen in both calculations to be a consequence of the vortex topology. A 
preliminary discussion is given regarding finite temperature extensions of the 
Berry phase calculation, with emphasis on its relevance for the sign anomaly 
occurring in Hall effect experiments on type-II superconductors in the flux 
flow regime. 

74.20.-z, 03.40.Gc, 74.60.Ge 



2 



Typeset using REVT^ 



I. INTRODUCTION 



As is well-known, the mixed state of a type-II superconductor is characterized by the 
partial penetration of magnetic flux into the superconductor in the form of flux lines, or 
vortices. Efforts to understand the dynamical behavior of such vortices have persisted for 
almost 30 years, and though much progress has been made, a number of basic issues continue 
to be controversial. A vortex in a type-II superconductor is acted on by 3 classes of forces: 
(i) pinning forces due to lattice defects; (ii) dissipative forces due to coupling of the vortex 
core to the lattice; and (iii) a non-dissipative force F^^ due to an applied magnetic field 
liext, an electric field E generated by the vortex motion, and the hydrodynamic pressure 
of the surrounding condensate of superconducting electrons. Already in the classic models 
of Bardeen-Stephen (BS)0 and Nozieres-Vinen (NV)i, the form of F^^ is controversial. The 
debate centers around whether F^^ includes a contribution from the Magnus (lift) force of 
classical hydrodynamics which acts on a solid body moving through a fluid which circu- 
lates about iti. Both models: (1) are macroscopic/phenomenological in character, based 
on Maxwell-London electrodynamics, classical thermodynamics, and physical intuition; (2) 
assume strongly type-II superconductors so that the non-local character of BCS supercon- 
ductivity can be approximated by a local dynamics; (3) assume T = so that normal 
electrons are only present inside the vortex core whose radius is equal to the zero tempera- 
ture coherence length ^o; and (4) the applied magnetic field satisfies Hc^ < Hext ^ so 
that vortex-vortex interactions can be ignored. The vortex is assumed to be immersed in 
an applied transport current J = Psev^, where ps is the superconducting electron density 
far from the vortex; e is the electron charge; and is the velocity of the applied super- 
current with respect to the lattice rest frame. In the BS model, the non-dissipative force 
(per unit length) is due to the Lorentz force, pshuoiys x z)/2. Here h is Planck's constant; 
and uj = ±1 is the vortex winding number whose sign specifies the sense of the condensate 
superflow about the vortex (which is threaded by a single flux quantum 0o = hc/2e). In 
this paper it will be assumed that all vortices are rectilinear (viz. line) vortices whose axis 
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lies along z. In the NV model, Fnd is the sum of the Lorentz force and the Magnus (lift) 
force —psTnlCvL x z, where m is the electron mass; /C = huj/2m is the circulation of the 
condensate near the vortex core; and is the translational velocity of the line vortex with 
respect to the lattice. The discrepancy in the form of Y^d in the two models leads to different 
predictions for the Hall angle in the flux flow regime^. Recent work by Ao and ThoulesJi 
(A&Th) has renewed interest in this controversy concerning F^^^. These authors argue that 
the correct form for Ynd is the NV-form, Ynd = Pshuj{ws — v^) x z/2, and that the Magnus 
(lift) force contribution arises from a Berry phase induced in the many-body ground state 
by the vortex motion. Their analysis is based on an ansatz for the many-body ground state 
and so it might be objected that their result concerning the presence of the Magnus force is a 
consequence of their ansatz and not an actual property of the true superconducting ground 
state. Thus it is of considerable interest to see if F^^^ can be determined on the basis of a 
microscopic analysis in which the superconducting dynamics is treated exactly. 

In this paper, two calculations of the non-dissipative force F^^ will be provided which are 
based on the microscopic formulation of the superconducting dynamics due to Bogoliubovi. 
This approach is powerful enough to treat problems with spatial inhomogeneities as occur 
when a vortex is present in the superconductor. In Section ^ F^^ will be determined by 
working with the true superconducting ground state in the case where a vortex is present. We 
first show how this state is constructed from the exact solutions of the Bogoliubov equations 
in the presence of a vortex. This state is seen to develop a Berry phase as a consequence 
of the vortex motion. Our approach in this subsection is a significant generalization of one 
used previously to determine the intrinsic orbital angular momentum in ^He-A for a spatially 
uniform orbital texturei. The Berry phase is calculated and found to agree with the result of 
Ref. ^ in the case of a neutral superfluid. A discussion of why the approximate calculation of 
A&Th has produced the exact result is also given. In the case of a charged superfluid, gauge 
invariance requires a correction to the Berry phase obtained in Ref. Q however this correction 
will be seen to not influence the final result for F^^. It is then shown how the Berry phase 
induces a Wess-Zumino term in the action describing the hydrodynamic degrees of freedom 



of the superconducting condensate. Variation of the hydrodynamic action with respect to 
the vortex trajectory gives F^^ which is found to have the NV-form. This calculation also 
shows that the Magnus force contribution to Fnd is a consequence of the Berry phase induced 
by the vortex motion. As our calculation is based on the true superconducting ground state, 
we see that the result of A&Th is not a consequence of their ansatz, but a real property 
of type-II superconductivity. In Section |T|, Fnd is determined via a microscopic derivation 
of the acceleration equation for the superflow. We find the expected contributions off the 
vortex due to spatial variation of the chemical potential, and to the electric and magnetic 
fields present. But by carefully tracking the effects of the vortex topology on the dynamics, 
we also find a singular term that describes the disappearance of linear momentum into the 
vortex, and which represents the non-dissipative force acting on the vortex. We find that this 
force is also given by the NV-form of F^^, in agreement with the Berry phase calculation 
of Section ||. As NV have shownll, the terms in the acceleration equation acting off the 
vortex lead to a flux of linear momentum in towards the vortex and corresponds to their 
determination of Fnd- By keeping track of topological effects, we will see that the rate 
at which linear momentum is appearing inside the vortex is exactly equal to the rate at 
which it is flowing in towards it so that linear momentum conservation is maintained in the 



combined condensate-vortex system as expected. We close in Section |^ with a discussion 
of our results; and with preliminary remarks concerning finite temperature effects and their 
possible consequences for the sign anomaly occurring in Hall effect experiments on type- 
II superconductors in the flux flow regime. Two Appendices are also included: Appendix 
A indicates how the Berry phase appears in a propogator formulation of the Adiabatic 
Theorem; while Appendix B gives a microscopic derivation of the gauge invariant Wess- 
Zumino term which is independent of the Berry phase calculation given in Section and 
also obtains the hydrodynamic action for the superconducting condensate. 
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II. BERRY'S PHASE AND THE NON-DISSIPATIVE FORCE Fnd 



This Section is organized as follows. In Section 11 A we construct the true supercon- 
ducting ground state in the case where a single vortex is present and derive an expression 
for the Berry phase induced in this state by the vortex motion in terms of the Berry phase 
induced in the eigenstates of the Bogoliubov equation (in the presence of a single vortex). 
In Section 11 B the Berry phase in the true superconducting ground state is evaluated. In 
Section 11 C we show how this Berry phase induces a topological Wess-Zumino term in the 
action describing the hydrodynamic degrees of freedom of the superconducting condensate. 
Variation of the hydrodynamic action with respect to the vortex trajectory gives F^^ and it 
is seen in the course of evaluating this variation that the Berry phase induced Wess-Zumino 
term is responsible for the occurence of the classical Magnus (lift) force in F^^. 



A. Constructing the Ground State 

Our starting point is the Bogoliubov equationil. This formulation of the quasiparticle 
dynamics is capable of handling the spatial inhomogeneities introduced by a line vortex 
which for the time being is assumed to be fixed at the origin (in the lattice rest frame). As 
in the models of BS and NV, we will: (i) assume T = so that the physical behavior of 
the superconductor is captured by the superconducting ground state; (ii) approximate the 
non-local character of the BCS superconducting ground state by a local dynamics which is, 
strictly speaking, only true for strongly type-II superconductors (it is not anticipated that 
non-local effects will qualitatively modify our results); and (iii) assume Hc^ < Hext ^ 
so that vortex-vortex interactions can be ignored and we can focus on a single vortex. We 
adopt cylindrical coordinates (r, 9, z) with the z-axis along the vortex axis. In the presence 
of a line vortex, the gap function takes the form A(r) = Ao(r) exp[— i^^] so that the phase of 
the gap 0(r) = —0{v) is not a single valued function of position. As the field point r encircles 
the vortex axis, (p ^ (p — 2tt (more generally, A(r) = Ao(r) exp[ico'6'] so that — + 2ttuj, 
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where uj is an integer known as the vortex winding number). The gap amphtude Ao(r) goes 
hnearly to zero at the origin and approaches a constant value far from the vortex. In the 
presence of a magnetic field H = V x A, the Bogoliubov equation is 

\ 

{En — Hbog) 
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Here Ef = p'j/2 is the Fermi energy {h 
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c = 1 unless explicitly stated otherwise); n 



is an index which labels the spectrum of Hbog] Uo{r) is the pinning potential which will 
henceforth be set equal to zero so that we are dealing with a clean type-II superconductor. As 
is well-known, equation (|I]) has positive and negative energy solutions (relative to the Fermi 
surface) which are related as follows: if («„ f„) is (the adjoint of) a solution of equation (|ip 
with energy En > 0, then (— t"* m*) is (the adjoint of) a solution of equation (|l|) with 
energy —En < 0. The two-component elementary excitations built up from the solutions 
of equation (|l]) will be referred to as Nambu quasiparticles (NQP) and quasi- holes (NQH). 
Let al l create a positive energy (_E„ > 0) NQP and 2 create a negative energy (_E„ < 0) 
NQP. Introducing the two-component field operator \E'(r) for a NQP, we expand it in terms 
of the solutions of eqn. (|l|) 

(3) 

\E'(r) can also be defined in terms of the field operators for the spin- up and spin-down Landau 
quasiparticles (LQP) 
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The right-most part of eqn. (|^) comes from expanding the LQP field operator in terms of the 
creation and annihalation operators { 7„|, 7^^; 7„j, 7^^ } for the Bogoliubov quasiparticles 



(BQP)i. Comparison of eqns. ® and @ indicates that a„,i = 7n| and a„,2 = 7„|- Thus the 
spin projection 2sz labels the charge conjugation degree of freedom of the NQP'si and the 
superconducting ground state is obtained by occupying the negative energy states 

isc5)=n<2io) = n7niio) , (5) 

n n 

where |0) is the zero-particle state. 

As Hbog is Hermitian, the solutions of eqn. (|1|) satisfy orthogonality and completeness 
relationsi. The orthogonality relations are: 

j d^r [Un{v)u*^{r) +Vn{v)v*^{v)] =5nm. 

d^[un{v)vrn{r) - Vn{r)um{r)] = ; (6) 



and the completeness relations are: 

E [n„(r)<(r') + <(rK(r')] = 5(r - r') 

n 

E [^n(r)<(r')-<(r)n„(r')] = . (7) 

n 

Use of the orthogonality relations allows us to invert eqn. yielding 

7ni = / d^'T [ -V^}(rX(r) + ^|(r)M*(r) 

-il^= I dh [ij\{v)Un{v)+i^^{v)Vn{v)] . (8) 



As we shall see in the following sub-section, adiabatic motion of the vortex produces a 
Berry phaseS in the solutions of eqn. (P, (m„ Vn) exp[z0„](u„ f„), where 0„ depends 
only on the vortex trajectory ro(t). From eqn. (|D we see that 7„| inherits this phase, 
Ini exp[— z0„]7„j. Consequently (see eqn. (|D), the superconducting ground state develops 
a Berry phase F, \BCS) — > exp[iF] ISCS), where 

r = -E<^" • (9) 

n 

The physical significance of F will be discussed below. First we evaluate {0„}. 
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B. Calculating the Berry Phase 

BerryH showed that when a quantum system is coupled to an adiabatically-evolving 
environment and the system is initially prepared in an eigenstate of the t = Hamiltonian 
i/[R(0)] (the adiabatic coupling of the system to its environment enters through the set 
of parameters R(t) = (-Ri(t), ■ ■ ■ , Rn{t)) which appear explicitly in the system Hamiltonian 
if[R(t)]), that the time development of this state acquires a non-integrable phase ■ysit) now 
known as the Berry phase 



\ij{t)) = exp 



ilBif) —if duE{u 
Jo 



\m) ■ (10) 



Here \E{t) ) is the eigenstate of the instantaneous Hamiltonian if [R(t)] which evolves con- 
tinuously from the initial state \E{0)), and the gauge invariant Berry phase appropriate for 
an electrically charged particle is given by0 

Mt) = £ driEl^±-lMr)lE)^ I' ,riE\.f.V^-lA„m . (11) 

In Appendix A we show how the Berry phase arises in a propogator approach to the Adi- 
abatic Theorem which will be useful for our discussion in Section [11 C| . In our case, the 
quantum system is the two-component NQP field the environment is the superconduct- 
ing condensate described by the gap function A(r); and the adiabatically- varying parameters 
correspond to the vortex position ro(t) = {xo{t), yoif) )■ To determine V in eqn. (P) we must 
obtain the eigenstates of Hbog for a given vortex position. These states were first considered 
by Caroli et. ahlli, and a detailed solution for the bound and scattering states was given by 
Bardeen et. al.i. These eigenstates take the form 



(u \ 

Here hkz is the z-component of the quasiparticle momentum; 2/i is required to be an odd- 
integer to insure single- valuedness of the eigenstate Xn] tan6'(r) = {y — yo)/{x — xq); and 



r = y {x — xqY + {y — VoY- The spinor /„(r) has been examined in Ref. ^ analytically for 
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a simple model of the gap function, and numerically via a variational method. We will not 
require its detailed form below. To simplify the notation, we have collected the quantum 
numbers specifying the eigenstate into the single label n. Finally, our calculation is done 
per unit length in the z-direction so that (Px — * (Px. Thus 

^.-^jd^.j dr ro ■ xi V.X^ ' j d^^ j dr^ Mr) xix. 
= J d^xdrto- Xn} Vro^ + /Hr)Vro/n(r) 

-J d^xdT^Ao{T)xixn 
= 4>l, + 4>l- j d^xdr^Ao{r)xixn • (13) 



Combining eqns. (|]) and (|T3|) gives 

r = -E {€ + <Pn)+ I d^xdr^-AoY^xlXn 

= -E + I d'xdr^AoPs 

n J a 

= r^ + r + j d^xdr^AoPs . (14) 

At this point it is necessary to recall a basic property of BCS superconductivity. Specifically, 
that the global symmetry corresponding to changing the phase of the LQP field operator is 
spontaneously broken. Because of this, the ground state becomes infinitely degenerate and is 
parameterized by a phase angle / (see ref. |T^. The phase / and particle number N become 
canonically conjugate observables and two representations are possible for the many-body 
states. In one representation, particle number is well-defined, but / is not; while in the 
second, / is well-defined, but particle number is not0. Our calculation corresponds to the 
latter case where the many-body states do not contain a definite number of particles. Because 
of this, a calculation of ground state properties will require what amounts to a normal 
ordering prescription. An example of the kind of normal ordering we have in mind is the 
fixing of the Fermi energy through the requirement that the expectation value of the particle 
number operator in the ground state equal the number of electrons present: J2k = 
N. Because our calculation of F is based on a ground state in which particle number is 
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not well-defined, a similar kind of normal ordering is necessary and can be implemented 
straightforwardly. 

We begin by re- writing as 

r = - J dTd^xro-{Vror)F{r) , (15) 

where F{r) = Y.n fn^fn- F{r) receives a divergent contribution from the scattering states 
present in the Fermi sea. This divergence is regulated by introducing a temporary cutoff 
in the sum over n which is removed at the end of the calculation. We now show that the 
finiteness of Freg{r) implies that = 0. Writing tq = Vr x C, where C = {xqu — yox)z, in 
eqn. (|l^) and using standard vector identities gives 

Keg = dr / dl Freg{r)\C\ fi ■ (z X e^) - / d'^xC- { Vr-Freg X Vr-r + FregVr X V^r) 
= , (16) 

(5*00 is the circle in the x-y plane at infinity) since: (i) n = e^; (ii) V rFregir) and VrT are 
parallel; and (iii) r is a single- valued function of (x, y). Removing the cutoff gives = 0. 
can be re-written as 



^' = -J drd^x {ro-VroO)Yl l«n|' (-/i + ^) + l^^nl' - ^) 



(17) 



= - I drd^x (i-o ■ Vro^) S . (18) 

From eqn. ( [TtI) and ([12|) , we see that 5* is minus the z-component of the orbital angular 
momentum density (in units of h) present in the Fermi sea. Although our expression for S 
is divergent because this sea is infinite, knowledge of its physical significance makes obvious 
what the physically relevant normal ordering choice is: 

Sreg = ( "A/^ ) = -Z ■ F X {psmv)/h . (19) 

Here v = ee(/C/27rr)/(r) is the superflow about the vortex; /C = huj/2m is the circulation 
of the condensate near the vortex core {uo is the vortex winding number, see above); and 
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1 Co < < A 
fir) = { . (20) 

~ re"''/^ r > A 

It is extremely important to realize that this form of /(r) insures that only one flux quantum 
threads the vortexlli. Thus our normal ordering choice, eqn. is fixed by the requirement 
that the flux through the vortex $ agree with our assumption $ = hc/2e. (In the following 
subsection the ground state Berry phase will be related to the occurrence of a Wess-Zumino 
term in the hydrodynamic action. For an independent derivation of this Wess-Zumino term 
which does not rely on the above normal ordering prescription, see Appendix B.) With this 
choice of normal ordering 

Sreg = -ps^/'^ , (21) 

and 

T = j dT(fxps{Uo-V,,d-^A,^ , (22) 

for C(j = — 1, as assumed in our calculation. Thus we see that the true superconducting 
ground state does in fact develop a Berry phase as a consequence of the vortex motion. 
We are able to obtain this Berry phase and we see that the approximate calculation of 
Ref. H has produced the exact result in the case of a neutral superfiuid. For a charged 
superfluid, the gauge invariant form of the Berry phase must be used which adds the "eAo"- 
term to the ground state Berry phase. It is clear from our calculation that, in the case 
of a neutral superfiuid, the Berry phases are ultimately due to the azimuthal symmetry of 
the vortex which requires the eigenstates of the Bogoliubov equation to be eigenfunctions 
of the z-component of orbital angular momentum. Examination of the ansatz of A&Th 
shows that it is such an eigenfunction and so has the essential ^-dependence necessary for 
a proper determination of the Berry phase in the case of a neutral superfiuid. Because we 
have worked with the true superconducting ground state it is clear that the gauge invariant 
Berry phase and its physical consequences, to be discussed below, are true properties of 
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the superconducting ground state when a moving vortex is present, and not the result of a 
specific ansatz wavefunction. 



C. Berry's Phase, the Hydrodynamic Action and F^^ 



Our starting point is the vacuum-to-vacuum transition amphtude for a system of elec- 
trons with an effective attractive interaction due to phonons which is responsible for the 
pairing instability which characterizes BCS superconductivity, 



W = {vac-t = T\i){t = T)) = {vac;t = T\U{T,0)\vac;t = 0) . 



(23) 



Here [/(T, 0) = T(exp[— i /q (ItHbcs])] ^ is the time ordering operator; and the BCS 
Hamiltonian is 



BCS 



d X Tpl^x 

_ 9_ 
2 



1 9 

-- (V - leAf ~Ef + eAo 



+ / d^x 



Stt 



fH — H„T. 



ext I 



(24) 



Here is the field operator for a Landau quasiparticle with spin a; H = V x A is the 

microscopic magnetic field and its associated vector potential; Hext is an applied magnetic 
field; E = —(VAq + dA/dt) is the electric field produced by the vortex motion; g > is the 
(attractive) coupling constant; and repeated indices are summed over. The gap function A 
is most conveniently introduced via a Hubbard-Stratonovitch transformation0 which allows 
eqn. (|^) to be written as a path integral over (A, A* ) in which the four Fermion interaction 
is replaced by the bilinear Fermion pairing interaction 



W 



V[A]V[A*] {vac]A{t = T)|?7a(T, 0)|mc; A(0) ) 



(25) 



Here Ua{T, 0) = T(exp[— i drHeff] ), where -f^e// = Hj + L,. + Lgm] and 



Hf 



1 9 

- (V - ieAY -Ef + eAo 



d^x tpl 

L ^ 

+ / d^x \A*ipiip^ + Aip\i;l 



^<7 



(26) 
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f ,3 1^1' 




/ d'^x - — - + 


[ d^x — 


/ 2g 


J 8tt L 



(H-He.0'-E2 ■ (27) 



Note that the vacuum state and -f^e// both depend on the particular "path" taken by A. As 
discussed earher, we are hmiting ourselves to "paths" whose time dependence corresponds 
to adiabatic motion of the line vortex (i. e. the vortex velocity satisfies |ro| ^ f/; see also 
Bardeen0). 

The action for the superconducting condensate S is given by 
= (mc; A(T)|f/A(T,0)|mc; A(0)) 

= {vac] A{T)\UA{n){tn+i,tn) ■ ■ ■ UA{k)itk+i,tk) ■ ■ ■ UA{o)(ti,to)\vac-, A{0)) , (28) 

where Sq is the action for the bulk degrees of freedom of the condensate; Shyd is the action 
for the hydrodynamic degrees of freedom; A{k) = A(tfc); = T; and to = 0. To make use 
of the Adiabatic Theorem it proves useful to define the instantaneous eigenstates \En{t) ) of 
the instantaneous effective Hamiltonian Heff{t). Appropriately inserting complete sets of 
instantaneous eigenstates gives 

g-<5o+5,,,) ^ (^ac; A(r)|f/A(„)(t„+i,t„) X 

^El^K^fc))(^^(i^fe)I^^A(fc)(4,4-l)|^n^(^fe-l))(^m(^fc-l)| [ X 
_ k=l l,m J 

UAio)ih,to)\vac;A{0)) . (29) 

From the Adiabatic Theorem, we know that f/A(o)(^i; ^o) evolves \vac; A(0) ) into an eigen- 
state of Heffiti) (see Appendix A) which will be denoted |£'o(ti); A(ti) ). Similarly, 
UA{k-i){tk,tk-i) evolves \Eo(tk-i); A{tk-i)) into the energy eigenstate of Heff{tk) denoted 
\Eo{tk); A{tk) ) for n>k> 2. Thus 



e 



-i'^'^+'^'^y'') = {vaC;A{T)\UAin){tn+l,tn)\Eo{tny,A{tn))x 

n 

n(^o(4); A(4)|f/A(fc)(4,4-i)|^o(tfc-i); A(tfc„i)) X 

k=l 

(Eo(ti);A(ti)|f/A(o)(ti,to)|t^«c;A(0)) . (30) 
14 



From the Adiabatic Theorem (see Appendix A), we have that (e = T/{n + 1)) 



Eo(tfc); A(fc)|f/A(fc)(4,tfc_i)|^o(tfc-i); A(fc - 1) ) 

= (Eo(tfc);A(fc)|e'^^^e-*^°(*'=)^|Eo(tfc);A(A;)) 



(31) 



Here 7^ is the time derivative of the gauge invariant generahzation of the Berry phase at 



time tk- The matrix element appearing in the final line of eqn. (|3T|) can be written as a path 
integral over the fermionic degrees of freedom which can be evaluated exactly since Hj is 
quadratic in the fermion field operators. This calculation is described in Appendix B and 
the result is given by exp[— iS'][], where 

lAP 



S = So — / d X 



2^7 



d xdr 



m 



^p,v2(x, t) + Ar(0)A; 



(32) 



Here Sq is the action for the bulk degrees of freedom, which is not of immediate interest to 
us, and so will not be written out explicitly; ps is the density of superconducting electrons 
at T = 0; Vs = — (/i/2m)[V0 + (2eA)/(/;,c)]; (p is the phase of the gap function; A^(0) is the 
electron density of states at the Fermi level; Aq = cAq + {h/2)dt4>; ^i-nd h, m and c have been 
re-instated. Thus from eqns. (^), (0) and (|32D, the hydrodynamic action is 



S, 



hyd 



dr 



-^7b + / d,^^ 



mps 



v.^ + A^(0)i^ + ^{(H-He.,)^-E2} 



(33) 



We see that the Berry phase induced in the superconducting ground state by the vortex 
motion has found its way into Shyd- Such contributions to low energy effective actions by 
Berry phases are well-known and the induced term is referred to as a Wess-Zumino termlli. 
Wess-Zumino terms are topological in origin and in our case it is the vortex topology which 
is responsible for the occurrence of such a term in Shyd- We are interested in a superflow 



^The result quoted in eqn. (32) does not include the contribution 5*1 found in Appendix B. This 
is because we will obtain Si in this subsection from the ground state Berry phase T (found in 
Section^). 
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which is the sum of an apphed steady superflow v = {h/2m)'V(3 and one that circulates 
with velocity Vcjrc(r; ro(^)) about the adiabatically moving vortex. The gap phase acts as 
the velocity potential for the superflow (p = j3 — 6{r] ro(t)), and 6 is the azimuthal angle 
about the vortex at ro(t). 

In the absence of pinning centers, the energy of the ground state will not depend on the 
location of the vortex so that the initial ground state \vac; A(0) ) will evolve adiabatically 
into the ground state of Heff(t) for all t of interest. Thus 75 is the ground state Berry phase 
r calculated in Section [II B| . Making use of eqn. (p2D gives 



hyd 



drd X 



+ 



mps,.2 



(34) 



, '-^rl + N{0)Al 

The terms in Shyd linear in Vro^' describe the coupling of the vortex to the applied superflow 
v; to the electric and magnetic field via {Aq, A); and to the superconducting condensate via 
the Berry phase F. Thus we re-write eqn. 



as 



S, 



hyd 



dr d^x 

Scond Sint 



cond 



—A + 2eAo^^ro - ro 
mc Ps , 



(35) 



The non-dissipative force acting on the vortex is given by the functional derivative of Sint 
with respect to yq. Carrying out the functional derivative gives 



^n. = -^(v-ro)xz + o[|^^ 



(36) 



for a vortex with winding number a; = — 1, which is the case we have been considering. The 
effect of the scalar and vector potential on the vortex is seen to contribute to higher order in 
where ^0 is the zero temperature coherence length and A is the London penetration 
depth. For arbitrary winding number, we have the general result 

Pshuj 



nd 



roj X z . 



(37) 
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This is the principal result for Section || and is seen to agree with the result of Ao and 
Thoulessi which was based on an ansatz for the many-body ground state. In this section 
we have worked with the true superconducting ground state (in the presence of a vortex) 
and found that the Berry phase generated in this ground state is responsible for producing 
a Magnus (lift) force on the vortex as argued by Ao and Thoulessi, and that Fnd is given 
by the NV-form. In the following section we will focus on the momentum flows occuring in 
the superconductor and from this analysis we will be able to determine the non-dissipative 



force on the vortex. The result obtained in Section |T| will be found to be identical to the 
result in eqn. (pTf). 

III. LINEAR MOMENTUM FLOW ANALYSIS AND Fnd 

In this section we will examine the flow of linear momentum in a superconductor when 
a moving line vortex is present. Our microscopic analysis will be based on the Bogoliubov 
equation. In the first subsection we will re-write the Bogoliubov dynamics in a pseudo- 
relativistic notation which proves convenient for calculational purposes. In Section [III B| 
we work out the continuity equation for the condensate linear momentum at T = 0. This 
equation is seen to contain a source term which is shown in Section [III C| to be non- vanishing 



and follows from the non-invariance of the measure in a path integral description of the 



NQP dynamics under a phase transformation of the NQP field operator. In Section [III C 
we evaluate the variation of the measure to obtain an explicit expression for the linear 
momentum source term. In Section P^IIDj we show that the continuity equation is simply the 



acceleration equation for the superfiow. The topological character of the vortex is seen to 
lead to a transfer of linear momentum from the condensate to the vortex which corresponds 
to the non-dissipative force Fnd- Our analysis shows that the result found for Fnd from 
tracking the momentum flows agrees exactly with the result of the Berry phase calculation. 
The contribution to Fnd rooted in the topology of the vortex and manifesting in the non- 
invariance of the measure is further seen to contribute the classical Magnus (lift) force to 



17 



A. Preliminaries 



Our starting point is the Bogoliubov equation given in eqn. ([1|) for the case where a 
vortex is present with winding number uj = —1 (viz. A = Ao(r) exp[—i6]). Substituting 



\"»/ 



62 



(38) 



into eqn. (|l|) leads to the Bogohubov Hamiltonian 
Hbog 



( \{iV-v,f-E 



\ 



0-3 



An 



/ ^0 



+ Ao (Ti 



(39) 



where the {ai} are the 2x2 Pauh matrices corresponding to the two-dimensional Nambu 
space and have nothing to do with spin; and Vs = —{1/2)'V6 — eA {h = m = c = 1). 
We now make an eikonal approximation for the eigenstates near the Fermi surface^ 



(u ^ 



(40) 



where |q| = kf and the (?/„, Vn) vary on a length scale kr^. To first order in gradients 



Hbog = (^3 [-q ■ («V - 0-3 vj] + Aq (Xi . 



This gives rise to the gauge invariant second quantized Lagrangian 



(41) 



idt + CTg ( -dtO - eAo ) + CTs q • {iV - (XgV,) - Aq CTi 



(42) 



We see that the original 3 + 1 NQP dynamics has separated into a collection of independent 
1 + 1 subsystems labeled by q = q/A;/ which will be referred to as q-channels. This separation 
allows us to reconstruct the 3 + 1 dimensional dynamics by adding together the contributions 
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from the separate q-channels. Because the dynamics of the q-channels is 1 + 1, and the Pauh 
matrices are the two-dimensional representation of the Dirac gamma matrices, we can re- 
write eqn. p^ ) in a pseudo-relativistic notation. With x° = t; = q ■ x; and 7° = 0"i, 



Y = -tCT2, 75 



(Ta; we can write 



(43) 



Here /i = 0, 1; = eAo — {l/2)dt9; Ai = q - v^; and ^ = \l/'^7°. Our final manipulation 
involves analytically continuing to imaginary time, x° = —ix'^, to obtain the Euclidean 
Lagrangian 



Here a = 1, 2; (Px 



-idPx = —idx'^dx^; Aq = A2, Ai = —iAi; and 7° = 7^, 7^ 



(44) 



27 , and 



75 = 75 = —27^7^. These definitions insure that the operator i'~)"'Da = —i'y"' {da + 75^0) is 
Hermitian. 



B. The Continuity Equation for the Condensate Linear Momentum 

The result of the last subsection is the separation of the 3 + 1 NQP dynamics into 
a collection of independent 1 + 1 dynamical subsystems to be referred to as q-channels. 
This arises from the eikonal approximation made for the eigenstates of Hbog near the 
Fermi surface in terms of wavepackets with mean momentum pj-q. By construction, both 
positive and negative energy eigenstates (viz. above and below the Fermi surface) carry mean 
momentum pj-q. 

Focusing on the q-channel, its single particle excitations will be quasiparticles with 
(mean) momentum pjq (right go-ers) and quasi-holes with (mean) momentum — P/q (left 
go-ers). In the superconducting phase, a Landau quasiparticle with (mean) momentum pj-q 
will be a superposition of right go-ers and left go-ers 

^q(a;) = e^^/^-^^K(x; q) + e-^^^^-''^i(x; q) , (45) 
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where 



V'i?(x; q)= J2 ^'^^'""(^k+ki 

k=-A 



(46) 



k=-A 



Qp (ap annihalates (creates) a LQP with momentum p; A is a cutoff that determines the 
spread of the wavepacket in momentum space; and spin indices have been suppressed. From 
this, the Nambu quasiparticle field operator is 



^q(x) 



^ ^r{^; q) ^ 



V 



^i(x;q 

The Noether current associated with the phase transformation exp[— is j^{x) 

^q7^\[^q. Using eqn. ( ^Tj) we see that (suppressing the q-dependence for the time being) 



(47) 



lim "^Ux + e)\&(x) 



lim 



V'|j(x + e)V'ij(x) + V^l(x + e)ipl{x) 



(48) 



Clearly, jo{x) is the net Fermion number density (particle minus hole density) and Pfq,jo{x) 
is the operator representation of the density of linear momentum along q 



[x;(i) =Pfqjo{x) 



(49) 



Taking the vacuum expectation value of eqn. (^91) and summing over all q-channels gives 
the density of linear momentum in the condensate 



9iW 



47r2 



cii{vac\f{x;ci)\vac)i 



(50) 



Here a is the spin index (±) and |fac)q is the vacuum (ground) state for the q-channel. 
Similarly, the condensate stress tensor is 



Tij = k)Y. J ^q^Qi (^«c|j^(x;q)|mc)q . 



(51) 
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The continuity equation for the condensate hnear momentum is then 

— cii{vac\dtf + cijdjj^\vac)ci 

= ^/E / ^qi(^«c|9^j'^|mc)q . (52) 

We go on to evaluate the matrix element appearing in eqn. and will see that it is non- 
vanishing due to the non-invariance of the measure in the path integral used to formulate 
the NQP dynamics. 



C. Non-invariance of the Measure 



The generating functional for the NQP Green's functions in the q-channel is 



(53) 



Here rj {rj) are external sources for \E' (\l/); and A/" is a normalization factor. Invariance of 
Wnqp under the change of variable \l/(x) — >• exp[— ?6'(x)]\E'(x) gives 

= j V[^] J d^y exp ij d'^x {i^d^ - Yl^A^ - Aq) ^ + r/^ + ^'r] } 

X e{y) -id^f + iF{y) +f]^ + ^rj . (54) 

F{ii) enters eqn. (Q) through the Jacobian of the transformation d[^'] = Jd[^] and (as 
we shall see shortly), J has the form: J = exp[z / d'^x 6{x)F{x)]. Setting rj = f] = in 
eqn. (Q) gives 

d^r = F{x) . (55) 
To calculate the Jacobian we follow Fujikawa0 and go over to Euclidean space x^ —ix^. 



From Section [III A| , this transformation leads to the Euclidean Lagrangian 



(56) 
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and 



iS 



(57) 



The continuation to Euclidean space has been done in such a manner as to insure that ija^a 
is Hermitian and so has a complete set of states i'y"'Da(f)n — ^n<l>n- 
Under the change of variable \1' exp [—i6{x)] ^f, 



where the field operator has been expanded in the complete set {0n}- Thus 



(58) 



(59) 



As is a Grassmann variable 



J = (det C) ^ = exp 
= exp 



d^XE (f)i{x)0{x)(f)n{x) 

n 

i J d^XEO{x)FE{x) 



(60) 



as anticipated, and 



Pe{x) = J2 (^i{x)(j)n{x) = J2 tr0n(a;)0n(^) ■ 



(61) 



Here tr is a sum over Nambu indices only. Fe{x) is divergent and its finite part is isolated 
by introducing a cutoff M and a subtraction to remove the divergent part. Thus 



Fe{x) = lim 

M2— >00 



(62) 



The {(f} are eigenfunctions of i'j°'Da\A=o = —ij'^da and {Xm} are the associated eigenvalues. 
Then, 



M2- 
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= -r^'^'Fab , (63) 

where = 1. In going from the second to the third hne the completeness of the {4>n} 
{{fm}) has been used as well as the 1 + 1 character of the q-channel dynamics. 
Switching back to Minkowski space gives F{x) = [e^'^ Fh^)/Ait so that 

d,f = ^ . (64) 

In the following subsection we will use this result to obtain the acceleration equation for the 
condensate superflow. 

D. The Acceleration Equation and F^^^ 



Inserting eq. ( plf ) into eqn. (m) , the continuity equation becomes 

= Co|-^[5o, 9,]0 + eE,| , (65) 

where Co = /cj/37r^ is the particle density in the normal phase when the chemical potential 
equals the Fermi energy; and h has been restored. The source term is clearly non-vanishing 
and contains the effects of the electric field whose origin is the vortex motion as well as a 
term whose origin is the topology of the vortex. The essential topological property of the 
vortex is that the phase of the gap function changes by 2tiuj as we wind once around the 
vortex (recall uj is the vortex winding number) so that 
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'c^x, dy](j) = 2'kuj5{x - XQ)5{y - yo) 



(66) 



where is the gap phase 



= —9 for our calculation). Thus 



We also note that since = —{h'V6)/2 — eA, 

h 



X (V X + -Vs X (V X V6') + ev^ X B = 



It is straightforward to show that (recall 



so that 



V, X (V X V^) = (v, X z) 5^{r - fq) 



Vs X (V X Vs) - ^ (v^ X z) 5^(r - ro) + ev^ X B = 



(67) 



(68) 



(69) 



(70) 



Eqn. (|70| ) is recognized as the (cross product of with the) London equation in the presence 
of a vortexH, and follows automatically from our definition of v^. Thus eqn. (|65D can be 
re-written as 



dtgi + djTij = C 



X (V X v^) + eE + ev^ X B - ^(v^ - tq) x z5^(r - tq) 



(71) 



It is possible to evaluate the LHS of eqn. ( [71] ) using eqns. (|5^) and (0) once the matrix 
element {vac\j'^\vac)q is known. A simple diagrammatic calculation along the lines of 
Goldstone and Wilczekli gives 



mc|j''|mc)q = —A^ . 

Zn 



(72) 



We will not reproduce that calculation here as eqn. ([7^ ) is clearly consistent with eqn. (|64D 
(though see ref. ^ if further details are desired). Thus 

d'^q ^ 1 ~ 



f f 
Co(v, 



(73) 
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and 



Co i^dtO - eAo] . (74) 



2 

Combining eqns. (O); O)) ^^^1 ([7^) ; together with some algebra gives 
d 

-Vs - X (V X V 



- V, X (V X V,) = (75) 



V(e/lo)-V(^9i^)+eE + ev, xB-:^(v,-ro)xz52(r-ro) . (76) 

From the Josephson equation {hdtcj)) = —fio, where /io is the chemical potential in the vortex 
rest frame and can be written as /io = /i + + cAq {fi is the chemical potential in the 
lattice rest frame and m = 1). Recalling that = —6, together with the preceeding remarks 
gives 

^v, + V(iv2) - V, X (V X V,) = -V/i + eE + ev, X B - ^ (v, - ro) x z ^^(r - fq) , 

(77) 

or finally, 

^ = -V/i + eE + ev, X B - ^ (v, - ro) x z ^^(r - ro) , (78) 

which is the acceleration equation for the condensate superflow. Thus we find the expected 
forces related to the hydrodynamic pressure (VP = p^V/i), and the electric and magnetic 
fields. We also see that linear momentum (recall m = 1) is disappearing from the condensate 
into the vortex at ro(t) at the rate {pshuj /2){\s — i"o) x z per unit length so that 

Pshoj 



nd 



ro)xz , (79) 



2 

in agreement with our Berry phase calculation. It is clear from eqns. (|66D, ( p7| ) and (p9| ) that 
the origin of this force is the vortex topology and that the Magnus force contribution to F„d 
arises from the non-invariance of the path integral measure. Our result is also consistent 
with Ref. 1^ which showed that the first 3 terms on the RHS of eqn. (|78D lead to a flux of 
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linear momentum in towards the vortex at the rate [pghuo /2){\s — i"o) x z which is exactly 
equal to the rate at which we find it appearing on the vortex. This topology driven fiow of 
linear momentum from the 3 + 1 dimensional condensate onto the 1 + 1 dimensional vortex is 
a particular example of the Callen-Harvey mechanism of anomaly cancellation0 (see below) . 

IV. DISCUSSION 

In this paper we have examined the non-dissipative force F^^^ acting on a line vortex in 
a type-II superconductor at T = 0. Our analysis throughout has been microscopic and is 
based on the dynamics inherent in the Bogoliubov equation. In the first half of this paper we 
have verified and strengthened an earlier result of Ao and Thoules^ which argued that F^^ 
contains a contribution from the Magnus (lift) force familiar from classical hydrodynamics 
and which they argued was a manifestation of a Berry phase induced in the many-body 
ground state by the vortex motion. We have been able to determine the Berry phase induced 
in the true superconducting ground state and find that the result of Ref. ^, which is based 
on an ansatz for the many-body ground state wavefunction, agrees with the exact result in 
the case of a neutral superfiuid. For a charged superconductor, gauge invariance requires 
an "ey4o"-term to appear in the ground state Berry phase. This correction is seen however 
to not effect the final result for F^^^. It is clear from our calculation that, in the case of a 
neutral superfiuid, the ground state Berry phase is ultimately a consequence of the topology 
and axial symmetry of the vortex. This requires the energy eigenfunctions to be eigenstates 
of the z-component of the angular momentum. Examination of the ansatz of Ref. |^ shows 
that it has captured this essential property of the true superconducting ground state, and as 
a consequence, their calculation produces the exact result in the case of a neutral superfiuid. 
Although our paper is clearly based on the BCS theory of superconductivity, the fact that 
the Berry phase is ultimately a consequence of the topology and axial symmetry of the 
vortex leads us to suspect that vortices in non-BCS type-II superconductors will likely 
induce Berry phases in their corresponding ground state. One would thus expect a similar 
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Magnus force contribution to Fnd in non-BCS type-II superconductors as well. By deriving 
the action for the hydrodynamic degrees of freedom of the superconducting condensate from 
the microscopic Bogoliubov dynamics we are able to see how the ground state Berry phase 
manifests itself as a Wess-Zumino term in this actionil. Variation of the hydrodynamic 
action with respect to the vortex trajectory gives the non-dissipative force acting on the 
vortex. This variation also shows that the Wess-Zumino term is responsible for producing 
the Magnus force contribution to Fnd- This calculation substantiates the discussion given 
by Ao and Thouless which again captures the essential issues involved. In the second part 
of this paper we have examined F^^ by examining microscopically the continuity equation 
for the condensate linear momentum. We show that this equation leads to the acceleration 
equation for the superflow and by being careful to track the effects of the vortex topology, are 
able to see that a flux of linear momentum flows from the 3 + 1 dimensional condensate onto 
the 1 + 1 dimensional vortex where it appears at a rate equal to {pshuo /2){ys — ro) x z (per 
unit length). The sink terms in the continuity equation for the condensate linear momentum 
are seen to arise from the topology of the vortex. This topology related transfer of quantum 
numbers (in our case, linear momentum) from a higher dimensional system Si in which a 
lower dimensional system 5*2 is embedded is an example of the Callen-Harvey mechanism 
of anomaly cancellation^. Specifically, linear momentum is not conserved in either the 
condensate or the vortex, but is conserved in the combined system. In our case, the issue 
of anomaly cancellation simply describes the linear momentum transfers expected between 
the condensate-vortex subsystems. It should be stressed that our calculational procedure 
is able to keep track of the more subtle momentum transfers that are a consequence of the 
vortex topology. Both calculations yield the same result for the non-dissipative force and 
verify the presence of the Magnus force in Ynd which is seen to be a manifestation of the 
vortex topology. 

It is straightforward to calculated the Berry phase for the many-body excited states. 
Positive energy NQP excitations are created with and positive energy NQH excitations 
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with 7^|. For a many-body excited state with NQP's in the single particle state (Pi, ■ ■ ■ , Pn) 
and NQH's in the states {Hi, ■ ■ ■ , Hm), the Berry phase is 

^MBiPi, ■■■,Pn; Hi,---, Hm) = - I] + (^^P^+ 

= r + re.e . (80) 

We see that a 2-fluid picture naturally develops in the Berry phase in which the contribution 
from the superconducting condensate is given by the ground state Berry phase V obtained in 
Section m together with a contribution Texc from the elementary excitations present. Because 
the Berry phase gives the orbital angular momentum along the vortex axis present in the 
energy eigenstates, we see that the total orbital angular momentum (along the vortex axis) in 
an excited state receives contributions from the condensate and the elementary excitations. 
From eqn. (^ we see that the orbital angular momentum associated with the elementary 
excitations is opposite to that of the condensate so that the orbital angular momentum 
of the excited states is less than the orbital angular momentum in the ground state (viz. 
Tmb < r). Because experiments are done at finite temperature, it is clearly of interest to 
generalize our T = analysis of F^^^ to finite temperature. An important question is the 
effect of temperature on the Berry phase contribution to the hydrodynamic free energy, and, 
consequently, of the effect of temperature on the Magnus force contribution to Fnd- Naively, 
one would expect this would lead to a thermal average of the Berry phase given in eqn. (pO|). 
As we have just seen, Tmb < T so that a thermal average of Tmb is expected to decrease 
with increasing temperature. Thus one would naively expect the Magnus force contribution 
to F„rf to correspondingly decrease with temperature. We should also mention that since 
r depends on the superconducting electron density which decreases and eventually goes 
to zero with increasing temperature, whereas Fg^cc increases with increasing temperature, 
it becomes clear that the thermal average of Tmb should change sign as a function of 
temperature. Because of the anticipated connection between the thermal average of Tmb 
and the Magnus force contribution to F^^, such a sign change will lead to a sign change 
in the Magnus force contribution to Fnd- Thus F^^ = {pshoj /2){ys — Oi{T)YQ) x z, and 
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(y(T) changes sign with temperature. Such a temperature dependence for a(T) will act, 
in general, to change the sign of the Hall angle in a Hall effect experiment on a type-II 
superconductor in the flux flow regime^. Such a sign change in the Hall angle as a function 
of temperature has been observed^ and the fact that neither the BS or NV models can 
explain this sign change is the reason for this effect being known as the "sign anomaly in 
the Hall effect" in superconductors. Further work on the temperature dependence of the 
Berry phase contribution to Fnd is in progress and is clearly of interest to efforts aimed at 
understanding the sign anomaly in the Hall effect. We will report on this work elsewhere. 
Volovikii has argued that, for temperatures satisfying A ^ fc^T ^ (A^/Ef), thermal 
broadening of single particle states bound to the vortex will lead to creation of excitations 
trapped on the vortex as a consequence of spectral flow and that this mechanism effectively 
wipes out the Magnus force contribution to F^^^ which has been found to be present at T = 0. 
Although Volovik's analysis is different from ours, it is clear that at these temperatures, the 
important excitations for a thermal average of Tmb will come from excitations bound to the 
vortex so that the remarks above regarding a weakening of the Magnus force contribution 
to Fnd are qualitatively in agreement with Volovik's result. It should be mentioned that 
Volovik does not discuss the possibility of the Magnus force contribution to F^^^ changing 
sign with temperature, nor of its possible consequences for Hall effect experiments on type-II 
superconductors in the flux flow regime. Ao et. al.0 have described a finite temperature 
calculation of the non-dissipative force and argue that the T = form of this force is 
preserved for T 7^ and that the only effect of temperature is to replace the zero temperature 
superconducting electron density with its finite temperature form. 

I would like to thank: Dr. Ping Ao for interesting me in this problem and for a number of 
helpful discussions; Dr. Michael Stone for interesting comments and discussions, particularly 
with regard to the issue of gauge invariance; T. Howell HI for constant support; and NSERC 
of Canada for financial support. 
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APPENDIX A: THE ADIABATIC THEOREM AND BERRY'S PHASE 

In this appendix we provide a derivation of the Adiabatic Theorem being careful not 
to overlook the appearance of the geometric/topological phase discovered by Berry0. Our 
analysis focuses on the appearance of the Berry phase in the propogator (viz. time devel- 
opment operator) for the wavefunctions as this is the relevant object for our analysis in 
Section y. For a standard treatment of the Adiabatic Theorem along these lines (minus the 
Berry phase), see Messiah^. 

As in our discussion of Section |T|, we consider a quantum system coupled to an envi- 
ronment which is evolving adiabatically. This coupling is described by the appearance of a 
set of parameters R(t) = (-Ri(t), ■ ■ ■ , i?Ar(t)) in the system Hamiltonian if[R(t)]. It proves 
convenient to introduce the instantaneous eigenstates of H[Il{t)], 



Clearly Ej{t) and \Ej{t) ) are continuous functions of R(t) so that one can introduce an 
operator A{t) such that 



For simplicity, assume the instantaneous energy eigenvalues are discrete and non-degenerate 
throughout < t < T. From eqn. (^), the projection operator Pj(t) = \Ej(t)){Ej(t)\ 
obeys 




(Al) 



E,{t))=Ait)\E,{0)) 



(A2) 



Pj{t) = AmmHt) 



(A3) 



Time evolution is governed by the propogator U(t,to) which satisfies 



zh-U{t,to)=H{t)U{t,to) 



(A4) 



We now state the Theorem to be proved: 



Adiabatic Theorem: 
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If initially the system is in the eigenstate \Ej{0) ) (of -ff[R(0)]), then the state of the system 
at time t is given by the eigenstate of if[R(t)] that evolves continuously from \Ej{0) ) (viz. 
\Ej{t) ) to within a phase factor); or formally, 



Ihn U{t,0)\E,{0))=Pj{t) \im U{t,0)\E^{0)) 



(A5) 



Proof: 



The operator A[t) is determined by the initial condition A[0) = 1 and the differential 
equation 



where K'{t) must satisfy 



zh— = K'{t)A{t) 



dP- 



(A6) 



(A7) 



Because of the presence of projection operators in eqn. (|A7|), this equation does not determine 
K'{t) uniquely. The most general solution can be shown to be 



K'it)=^hY: 



+Pf 



P . 



(A8) 



Here {fi{t)} is an arbitrary collection of operators. Requiring that the dynamics satisfy the 
Schrodinger equation will allow us to fix the {fiif)}. 

In the adiabatic limit f/(t, 0) A(t)(^{t) as t — ^ oo, where 



(A9) 



In general, f/(t, 0) = A{t)^{t)W{t), where W(t) is determined by the condition W{Q) = 1 
together with 

dW i 



dt h 



K'W 



(AlO) 



and 



K' = ^^A^K'A^ . 
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(All) 



By sandwiching eqn. (|A10|) between the states {Ej{0)\ and \Ei{0)), it can be shown that 
Wij = {Ej{0)\W\Ei{0)) ~ C(l/T) {i ^ j) and so vanishes in the adiabatic hmit. The 
diagonal elements are easily shown to be given by 



Wu = { E,{0)\W\Ei{0) ) = e^/o "^^^^ 



Thus 



i 

Thus for the situation envisioned in the premise of the Adiabatic Theorem, 

\m) = uit,omio)) 

This proves the Adiabatic Theorem. To fix the phase, we now determine K-^ 
By definition. 



It 
1 



K'\E,{t)) 



\E,it)) 



dt 



+PofA\m) 



Taking the time derivative of the definition Pj{t)\Ej{t) ) = \Ej{t) ) gives 



d 



d 



Equating eqns. ( [A15| ) and ( |A16| ) gives 



P,J)\E,{t)) = P,-\E,{t)) . 



(A12) 



(A13) 



(A14) 



(A15) 



(A16) 



(A17) 
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From this we find 



E,{t)\f,\E,{t))=z{E,{t)\-\E,{t)) 



and so 



'E,{t)\j^\m)+Mm\^\m) 



But the second term on the RHS of eqn. ( |A19|) vanishes since ^ [ ( Ej{t)\Ej(t) ) 
K'^^=th{EMi\EAt))^^d 



(A18) 

(A19) 
0. Thus 



(A20) 



Thus we recover the Berry phase contribution to the phase factor appearing in the Adiabatic 
Theorem. 

The derivation of the Berry phase given above assumes an electrically neutral system. 
For a system composed of electrically charged particles such as the electron gas in a su- 
perconductor, the Berry phase must be generalized in a manner consistent with the gauge 
symmetry underlying electromagnetism. This is because the hydrodynamic action Shyd must 
be gauge invariant, and the Berry phase induced in the superconducting ground state finds 
its way into Shyd (as seen in Section II C ). The appropriate generalization has been given 
by Aharanov and Anadan0. We will simply quote their result as the correct form can be 
guessed straightforwardly from knowledge of how the Berry phase transforms under a time 
dependent gauge transformation. 



' dr{E{T)\^4z-TMr)\E{r)) 



dr h 



(A21) 



Eqn. (|A2 1|) will be used in Section || of this paper. 



APPENDIX B: THE HYDRODYNAMIC ACTION REVISITED 



In this Appendix an alternative derivation of the gauge invariant Wess-Zumino term 
will be given. The calculation to be described is based on Eckern et. al.@. We will briefly 
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present their approach and show how the Wess-Zumino term arises very naturally in this 
approach. The result found will be seen to agree with the Berry phase calculation presented 
in Section [11 C| . It is hoped that the independent calculation of the Wess-Zumino term given 
in this Appendix will reassure the reader that no sleight-of-hand was perpetrated in the 
normal ordering of the Berry phase in Section [11 B| . A secondary aim of this Appendix is 
to point out that the effective action for a single superconductor obtained in Ref. |18| also 
contains a gauge invariant Wess-Zumino term (viz. a gauge invariant term first order in time 
derivatives of the gap phase). 

As in Section [11 Q , our starting point is the vacuum-to-vacuum transition amplitude 
{h = m = c = 1) 

rT 



W 



\ vac; t = T\T I exp 



-t I dtHBcsit) 

10 



\vac\ t = 0) 



(Bl) 



and 



BCS 



d^xtpl^x) 



|y d^xi)l{x)ipl^{x)ip^^{x)i)^{x) 



-- (V - ieAf -Ef + eAo 



d^X 



1 



57r 



fH-H 



ext I 



(B2) 



(See Section [11 C[ for a definition of the symbols.) The quartic interaction term is removed 
via a Hubbard-Stratonovitch transformationB so that W is expressed as a path integral over 
(A, A*) 



W 
where 



y" I?[A] I?[A*] (mc;t = T| T l^exp -i£dtH^ff{t) j 



\vac] t = 0) 



(B3) 



eff 



e 



(B4) 



and 



Hf 



d^xtpl{x) 



1 9 

- (V - icAy -Ef + eAo 



+ j d^x \A{x)i)\{x)i)\{x) + /\*{x)i)i{x)i)i^{x) 



(B5) 
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,3 W r ,3 1 

d X - — - + / d X — 

2g J 8tt 



fH-H 



ext ) 



(B6) 



The condensate effective action S is defined through the matrix element appearing in 
eqn. (pi) 



dtHf{t) 



|fac;t = 0) . (B7) 



The matrix element in eqn. ( |B7| ) can be calculated easily once Hf is re-written in terms of 
the NQP field operator ^ introduced in Section ||, 



where 



Hf = j d^xd^x''^\x)HBOG{x-x')^{x') ■ 



Hbog{x] x') = { (V - lea-^Af - Ej + eA^ | trg + Re{/\)a^ - Im{A)a2 



(B8) 



2 

X 6^^ - x')6{t - t') 



(B9) 



With Hf given by eqn. ( p8| 
( vac; t = T\T exp 



-I j dtHf{t) 



\vac; t = 0) = exp 



Tr InG 



-1 



(BIO) 



where 



= tS'i^ - ^')dt5it - t') - Hbog 



(Bll) 



As in Section pi A| , we make ImA = through the unitary transformation Hbog 
exp[z(0/2)o-3] Hbog exp[-i{(f)/2)as]. This gives 



d_ /I. 



^i: + ( + ) ^3 - Aofxi 



5=^(x-x')5(t-t') 



+ 



Aqo-s + zv, ■ V - 2^sO"3 



5=^(x-x')5(t-t') . 



(B12) 



Here An and v, are defined as in Section [111 A 



and 5G-^ = ^G^/ + 6G;^^ + dG'^^ 



contains 3 contributions which are gauge invariant and are (respectively) first order in time 
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derivatives of the gap phase, first order in space derivatives (of the gap phase), and second 
order in space derivatives. Putting these results together gives 



S = iTi InG + Sc + Sem , 
{Sc and Sem are the time integrals of Lc and Lem respectively). 



(B13) 



The trace in eqn. ( B13 ) is evaluated perturbatively 



Tr InG 



-1 



Tr In ( G 



Y-l 



-1 



6G 



Tr InGo ^ + TtGqSG-' - -Ti Go6G-'Go6G-' + 

2 



(B14) 



For the purpose of finding the hydrodynamic action, it is only necessary to work out the 
expansion in eqn. ( pi4|) to second order in space/time derivatives of the gap phase. Thus 



S = So + S^ + S2 + ■■■ 
=So+Shyd+■■■ . (B15) 

Eckern et. aklll have determined 5*0 and 5*2. 5*0 is the action for the bulk degrees of freedom. 
It is not of immediate interest to us, and so will not be discussed further. 5*2 is found to be 

m 



So 



dt drx 



(B16) 



where m has been re-instated; ps is the density of superconducting electrons at T = 0; and 
N{0) is the electron density of states at the Fermi surface. 

It is our purpose to show that Si is the gauge invariant Wess-Zumino term found in 
Section ffTO, where 



S^ = Sa + Ssi=i Tr Go (^G^/ + 6G;, 



(B17) 



We will evaluate Sa below. A similar calculation (which will not be reproduced here) shows 
that Ssi vanishes. Thus we focus on 



Sti = Tr Go6G^ 



-1 
ti 



i J £xd^ytT {x\Go\y){y\5Gt^^\x) 

i j d^x d^y tr {Go{x - y) [io(i/)(T3 ^^(x - y)5(t^ - ty)] } 



(B18) 
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and tr is a sum over Nambu indices only. We are interested in evaluating eqn. ( |B18| ) to 
lowest order in gradients of the gap phase. As Aq is already first order in gradients, we can 
write 



io(y) = io(R-p/2)^io(R) 



where R = (x + y)/2; and p = x — y. Also, Go{x — y) can be written as 



(B19) 



Gq{x - y) 



Jk{x-y) 



(B20) 



Due to the vortex, Aq varies with position near the vortex core. In the local limit, Aq is 
constant except on a set of measure zero (viz. on the vortex core). We will approximate Aq 
by its constant value away from the vortex core. Putting all this together gives 



Sti = I d-'Rd-'pitr 
- I d^RAniR 



d^k 



^ik{x-y) 



(27r)4 A;o-((kV2)-E^)a3-Aoai 
^ tr[-zG'o(k)a3] 



io(R) as 



d^k' 



-ik' (x—y) 



(271)^ 



h 



dtO - eAo 



(B21) 



where U has been re-instated. Thus, 5*41 is the gauge invariant Wess-Zumino term found in 
Section P C|. Since Ssi = 0, we have 5*1 = Sa and 



hyd 



d^x dt Ps 



\dte - eAo 



+ / d^xdt 



2 Ps^i + N{0) -dtO - eAo 



(B22) 



We stress that the calculation of the Wess-Zumino term given in this Appendix is in- 



dependent of the Berry phase calculation given in Section |II Q and thus acts to confirm its 
result, which is that Shyd contains a gauge invariant term first order in time derivatives of 



the gap phase. The occurrence of this term in S appears to have been overlooked in Ref. [18. 
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